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I .  INTRODUCTION 


The  present  paper  is  concerned  with  the  theoretical  description  of  steady- 
state  deflagration  of  uncatalyzed,  homogeneous  solid  propellants.  Our  principal 
aim  is  to  seek  a  minimum  level  description  of  the  burning  process  which  can 
elucidate  the  relationship  between  the  fundamental  chemical  and  physical  pro¬ 
perties  of  a  propellant  and  such  combustion  characteristics  as  the  pressure 
and  temperature  sensitivities.  It  is  our  feeling  that  such  a  description  must 
involve  processes  occuring  in  both  the  condensed  and  gaseous  phases  and  account 
for  their  mutual  interplay  at  the  phase  boundary.  Accordingly,  the  basic 
idealization  of  the  burning  process  to  be  considered  in  this  paper  is  that  of 
a  single  gas-phase  chemical  reaction  and  a  single  condensed-phase  process 
responsible  for  surface  gasification.  Both  sublimation  (or  evaporation  from  a 
melt  layer)  and  pyrolysis  into  gaseous  products  are  considered  as  candidate 
possibilities  for  this  gasification  process. 

The  chosen  idealization  is  by  no  means  novel;  however,  the  considerable 
literature  on  the  subject  is  marked  by  diversity  in  notation,  approach,  and 
intent.  Comparison  of  burning  rate  predictions  is  hampered  by  these 
characteristics  and,  more  importantly,  by  the  fact  that  many  theories  are  not 
strictly  expressible  in  closed  form.  In  this  paper  we  derive  approximate 
closed  form  expressions  for  the  burning  rate  pressure  and  temperature 
dependences  based  on  what  we  perceive  to  be  the  salient  assumptions  of  a 
number  of  these  theories.  Within  the  mathematical  framework  adopted  for  this 
study  most  theories  fall  into  one  of  three  broad  categories  of  heat  feedback 
approximation.  Naturally,  some  theories  fit  into  our  scheme  better  than 
others,  hence  no  judgement  is  implied  regarding  those  theories  lightly 
touched  on  or  not  mentioned  at  all. 

Progress  toward  the  understanding  of  complex  phenomena  has  often  been 
made  through  the  use  of  idealizations.  It  is  unfortunate  that  the  subject 
idealization,  primitive  as  it  may  seem,  appears  not  to  be  amenable  to  exact 
solution.  Expressing  the  gas-phase  reaction  rate  in  an  Arrhenius  form  (as 
is  usual)  introduces  extreme  non-linearities  into  the  conservation  equations, 
making  even  numerical  solutions  far  from  trivial.  Notwithstanding  this 
fact  and  the  enormous  potential  detail  which  might  be  required  to 
characterize  the  gas-phase  processes  for  each  propellant,  the  observed 
burning  rate  as  a  function  of  both  pressure  and  initial  temperature  displays 
a  remarkably  simple  and  universal  behavior.  This  suggests  that  insight 
might  be  gained  regarding  at  least  the  broad  features  of  the  combustion 
phenomenology  through  development  of  an  appropriate  idealization.  The 
present  study  represents  an  attempt  to  determine  the  adequacy  of  existing 
notions  for  this  purpose. 


II.  FORMULATION  OF  THE  1-DIMENSIONAL  STEADY-STATE  PROBLEM 
A.  General  Approach 

At  our  chosen  level  of  idealization,  the  solid  propellant  burning  process 
consists  of  one  mole  of  unburnt  propellant  molecules  A  in  the  solid  undergoing 
a  physical  or  chemical  transformation  into  N^  moles  of  gaseous  B  molecules. 
Each  mole  of  B  molecules  then  reacts  chemically  in  the  gas  phase  to  produce  N2 
moles  of  product  molecules  C.  The  gas-phase  reactant  B  may  be  either  the  same 
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chemical  species  as  A  (differing  only  by  a  phase  change]  or  a  different  chemical 
species  than  A.  Figure  1  represents  these  processes  schematically.  In  principle, 
the  problem  can  be  approached  by  solving  the  steady-state  conservation  equations 
in  the  solid  and  gas  phases  separately  and  then  matching  the  solutions  at  the 
gas-solid  interface  in  a  manner  consistent  with  the  assumed  solid-to-gas  decom¬ 
position  mechanism. 

B.  Gas  Phase  Description 

The  gas-phase  problem  consists  of  only  two  different  chemical  species, 
denoted  by  B  and  C.  It  is  supposed  that  B  undergoes  reactive  conversion  to  C 
with  the  consequent  release  of  heat.  Our  problem  is  to  compute  the  conductive 
heat  feedback  to  the  surface  of  the  solid  under  constant  pressure  conditions. 

The  radiative  component  of  the  heat  feedback  is  assumed  to  be  negligible.  As 
mentioned,  if  one  expresses  the  reaction  rate  in  an  Arrhenius  form,  the  con¬ 
servation  equations  take  on  a  highly  non-linear  character,  placing  an  exact 
solution  beyond  reach.  For  this  reason  some  investigators  have  obtained 
approximate  solutions  in  an  attempt  to  gain  an  appreciation  for  the  inter¬ 
relation  of  physical  and  chemical  processes.  However,  even  these  approximate 
solutions  are  sometimes  so  mathematically  involved  as  to  endanger  their  purpose. 
Such  analytic  intricacy  inhibits  comparison  of  the  salient  features  of  these 
theories  as  well.  It  is  our  intent  to  re-examine  this  approach  by  expressing 
the  essential  elements  of  a  number  of  these  past  approximations  in  a  single 
framework  which  hopefully  clarifies  their  physical  implications. 

The  conservation  equations  for  the  steady-stap  include  those  for  mass, 
chemical  species,  and  energy  given  respectively  as1 * * 


d_ 

dx 


(pu)  =  0 


Mass  Conservation 


d_ 

dx 


(mBpVB)  +  R(x)  =  0 


Chemical  Species  Conservation 


31  »£>  -  ST  (VhBVB  *  VW  '  0  Energy  Conservation 


where  p  is  the  local  mass  density  of  the  gas  mixture,  u  is  the  local  velocity 
of  mass  motion,  mB  and  mc  are  the  local  mass  fractions  of  species  B  and  C,  VB 
and  Vc  are  the  local  species  velocities  (combined  mass  motion  and  diffusional 
components),  R(x)  is  the  mass  of  B  undergoing  reaction  per  unit  volume  per  unit 
time  at  x,  X  is  the  local  thermal  conductivity  of  the  mixture,  T  is  the  local 
temperature,  and  hB>c  is  the  local  enthalpy  per  unit  mass  for  species  B  or  C. 

A  chemical  species  equation  can  be  written  analogously  for  me-  The  mass  con- 


1 

R.D.  Geckler,  "The  Mechanism  of  Combustion  of  Solid  Propellants Selected 

Combustion  Problems ,  AGARD,  Butterworths  (London),  1954.  p.  289. 
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Figure  1 


Schematic  representation  of  processes  considered  in  burning  of 
idealized  solid  propellant. 
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servation  equation  can  be  integrated  to  obtain 


pu  =  constant  =  M  (1) 

where  M  is  sometimes  called  the  mass  regression  rate  and  is  related  to  the 
linear  regression  rate  r  through  the  solid  density  by  the  equation 


M  =  p  r. 

s  * 


(la) 


If  we  assume  negligible  particle  diffusion  by  thermal  gradient  (Soret 
effect),  the  species  mass  fluxes  can  be  written 


dmB 

mBpVB  =  mBpU  -  Dp  dF 


(2) 


and  analogously  for  m^,  where  D  is  the  diffusion  coefficient  (Dg  =  Dc  =  D)  . 
The  species  conservation  equation  can  then  be  written 


j 

F  (Dp  dF>  -  M  dF  -  "  0 


(3) 


and  the  energy  conservation  equation  becomes 


d  ,,dT.  „„  dT  „  -  dmB  dT  ,  C  _B. 

dx  (xdiJ  -  MCP  d?  *  R  =  pD  zr  ar  tcP  -  Cp> 


where  the  specific  heat  of  the  mixture,  Cp,  is  related  to  the  component  specific 
heats  CB  and  Cp  by  Cp  =  mBcj*  +  mcCp.  We  now  make  the  first  substantive  assump¬ 
tions, 

* 

d  r 

C  =  C  =  C  =  constant 
P  P  P 

X  =  constant 

With  the  first  of  these  assumptions,  the  energy  equation  becomes 


(x4h  -  MC  +  q(x)  =  0 
dx  v  dx^  p  dx 


10 


(4) 


where 


qW  =  (hB  -  hc)  R(x)  =  qgR(x) 


(5) 


o  o 

the  symbol  Qg  =  hg  -  hQ  being  used  to  represent  the  heat  of  reaction  per  unit 
reactant  mass  in  the  gas  phase  at  the  reference  thermodynamic  state  (denoted 
by  the  zero  superscripts) .  Making  use  of  the  X  =  constant  assumption  and 
imposing  the  boundary  condition 


dT 

dx| 


X=" 


0 


Eqn.  (4)  may  be  integrated  by  the  integrating  factor  method  to  give  the  con¬ 
ductive  heat  feedback  flux  <J> 

b 


5  xdx,  =  f  exP x)  qOOdx.  (6) 

A  ''O 

1  x=0 


As  it  stands  Eqn.  (6)  is  purely  a  formal  result  since  q(x)  is  only  known 
after  the  problem  has  been  solved;  however,  it  is  a  particularly  meaningful 
equation  from  a  physical  standpoint.  The  exponential  factor 


MGD 

exP(-  y  X) 


is  a  measure  of  the  relative  effectiveness  with  which  heat  released  at  some  dis¬ 
tance  x  from  the  surface  contributes  to  the  overall  heat  feedback.  As  might  be 
expected,  the  larger  X  is,  the  more  effective  is  the  heat  conduction  back  to  the 
surface.  The  larger  Cp  is,  the  more  heat  is  retained  at  x  and  the  less  heat  is 
returned  to  the  surface.  The  larger  M  is,  the  more  difficult  it  is  for  heat  to 
return  to  the  surface  against  the  convective  counter-flow. 

The  gas-phase  heat  release  function  q(x)  can  be  related  to  the  physical 
and  kinetic  parameters  of  the  gas  phase  as  follows: 


q (x)  =  QgR(x)  =  qg  (mgp) v  Rg 


(7) 


where  v  is  the  reaction  order  and  kc  is  the  reaction  rate  constant  assumed  to 
have  the  form  (R  being  the  universal  gas  constant) 


kG  = 


AGe-EG/RJ 


(8) 


11 


Expressing  p  in  terms  of  the  ideal  gas  law. 


P 


WP 

RT 


(9) 


where  W  is  the  mean  molecular  weight  of  the  gas  mixture.  As  previously  noted, 
N2  moles  of  C  appear  for  every  mole  of  B  that  reacts  so  that  the  molecular 
weight  of  C,  Wq,  is  related  to  the  molecular  weight  of  B,  Wg,  by  Wq  =  Wg/N2- 
The  number  of  moles  of  B  and  C  at  any  x,  ng  and  n^,  is  given  by 


mBP 


nB  Wr 


mcp  ^  (l-ihB)pN2 


The  average  molecular  weight  at  this  point  in  the  gas  can  therefore  be  shown 
to  be 


W„ 


W  = 


N2-mB(N2-l) 


(10) 


Combining  Eqns.  (7),  (8),  (9),  and  (10), 


q(x)  = 


%  j  K  [ 'X 

RV  TV  [N2-mB(N2-l)]V 


(11) 


where  it  is  understood  that  mB  and  T  are  functions  of  x. 

A  number  of  other  useful  relationships  can  be  obtained  by  making  the 
further  assumption  that  the  Lewis  number  is  unity,  in  which  case 


Physically,  this  assumption  links  the  heat  transfer  by  conduction  to  that 
transported  by  diffusion.  Mathematically,  Eqns.  (3)  and  (4)  are  no  longer 
independent,  i.e. ,  the  solution  of  one  can  be  related  to  the  solution  of  the 
other.  Using  Eqn.  (12)  one  can  show2  that  the  enthalpy  is  uniform  throughout 
the  gas  phase,  i.e.,  if  we  follow  a  unit  mass  sample  from  the  surface  at 
temperature  Ts  to  the  final  flame  temperature  Tf,  the  enthalpy  of  the  sample 


2D.B.  Spalding ,  "The  Theory  of  Burning  of  Solid  and  Liquid  Propellants," 
Combustion  &  Flame,  Vol.  4,  p.  59  (1960). 
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at  any  point  is  given  by 


mB (T)  hB(T)  +  mc(T)  hc(T)  =  hc(Tf) 


(13) 


where  hg  is  the  enthalpy  of  B  per  unit  mass  of  B,  etc.  By  our  previous 
assumption  that  Cp  =  =  Cp  =  constant,  we  can  write,  for  example. 


hB  "  hB  *  Cp  (T-T°>- 


(14) 


Substituting  Eqn.  (14)  and  the  analogous  expression  for  into  Eqn.  (13)  and 

making  use  of  the  constraint  m  +  m  =  1,  one  can  show  that 

B  L 


mB(T) 


Cp(Tf-T) 


(15) 


which  relates  the  solution  of  Eqn.  (3)  to  that  of  Eqn.  (4)  at  all  x  in  the  gas 
phase . 

C .  Solid  Phase  Description 


In  this  section  the  external  energy  flux  required  by  the  solid  in  order 
to  maintain  a  given  steady-state  M  and  Ts  will  be  determined.  We  shall  neglect 
any  radiative  energy  loss  from  the  solid.  (See  Section  IVC  for  futher  comment 
on  this  assumption).  It  is  envisioned  that  the  solid-to-gas  decomposition 
comprises  the  following  sequence  of  steps.  At  x=-°°  the  solid  is  entirely  A 
so  that  m^=l.  In  the  case  of  a  pyrolysis  reaction  it  is  assumed  that  a  unit 
mass  of  A  will  be  converted  into  mg”0  unit  mass  of  B  and  (l-mg~°)  unit  mass  of 
C.  Thus  at  x=-o,  i.e.,  at  the  condensed-phase  side  of  the  interface,  m^=0,  mg= 
mg~°,  and  mc=(l-mg’°).  At  this  stage  B  and  C  are  considered  to  be  adsorbed 
onto  the  surface.  Since  the  symbols  hg  and  h^  refer  to  the  gaseous  phase  of  B 
and  C,  we  shall  denote  the  enthalpy  of  B  and  C  in  the  absorbed  state  by  hg  s 
and  h q  s.  Taking  mg~°  /  1  in  general  allows  us  to  consider  the  case  where 
both  reactive  as  well  as  non-reactive  products  are  formed  in  the  condensed- 
phase  pyrolysis  reaction.  For  simplicity  we  have  assumed  that  these  non-reactive 
products  have  the  same  properties  as  C.  B  is  assumed  to  react  only  in  the  gas 
phase. 

The  steady-state  energy  conservation  equation  in  the  condensed  phase, 
neglecting  particle  diffusion,  may  be  written 

k  (s  s)  -  MCS  E  *  <»,<*>  -  0  <16> 

where  \  is  the  thermal  conductivity  of  the  solid,  Cs  its  specific  heat,  and 
qs(x)  t^e  heat  absorbed  by  the  solid  per  unit  volume  per  unit  time  from 
reactions  in  the  solid.  Cs  is  assumed  to  be  constant.  Integrating  Eqn.  (16) 
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from  x  = 


_oo 


to  x=-0  with  the  boundary  condition 


=0,  we  obtain 


X=-oo 


♦s'  5  S|  =  MCs<VTo>  -  M  lhl  -  V” 

1  X=-0 


h° 

B,s 


(1-m. 


-o 


(17) 


The  last  term  in  Eqn.  (17)  may  be  obtained  by  constructing  q  (x)  by  analogy  to 
Eqn.  (5)  and  integrating  the  species  conservation  equation  o§  m^  (similar  to 
Eqn.  (3)  but  with  zero  diffusion)  from  x=-°°  to  x=-0. 

As  noted  by  Buckmaster,  et  al.,^  there  is  an  additonal  energy  flux  require¬ 
ment  arising  from  the  desorption  of  B  and  C  at  the  interface.  In  our  notation 
it  is  given  by 

♦i  ■  mi°  iW-YstVJ  +  M  “-O  thc(T»)-hc,stTs)1  • 

Since  hg  S(T)  =  hg  s  +  CS(T-T°)  and  hg(T)  =  hg  +  Cp(T-T°),  etc.,  the  total 
energy  fiux  requirement  can  be  shown  to  be 

♦s  ■  ♦s’  +  *i  ■  MW  r-To  -  h  <18> 

p  p 

where  Qs  =  -  mB~°h§  -  (l-mg"°)h®  +  (Cs-Cp)T°.  <j>s  in  Eqn.  (18)  simply  repre¬ 

sents  the  heat  flux  required  from  sources  external  to  the  solid  to  maintain  a 
surface  temperature  Ts  when  the  mass  regression  rate  is  M  and  the  initial 
temperature  is  T0.  <f)Q,  on  the  other  hand,  is  the  heat  flux  which  the  gas- 
phase  processes  are  capable  of  providing  at  a  given  pressure  when  the  mass 
flux  is  M  and  the  surface  temperature  is  Ts. 

We  shall  also  wish  to  examine  a  simple  phase  change  process  as  the 
mechanism  for  surface  gasification.  In  this  case  mg"°  =  1  and  Eqn.  (18)  is 
changed  only  by  the  substitution  Q  =-Q  ,  where  Qv  is  the  latent  heat  of  trans¬ 
formation. 

D.  Surface  Decomposition  Mechanism 

In  addition  to  having  to  solve  the  conservation  equations  for  the  gas  and 
solid  phases  separately,  one  must  also  have  a  quantitative  description  of  the 
mechanism  by  which  the  solid  decomposes  into  a  gas  at  the  interface.  The  de¬ 
tails  of  this  process  are  poorly  understood  despite  its  central  importance  to 
propellant  combustion.  There  are,  however,  two  schools  of  thought  which  have 
found  recurring  application.  One  is  that  the  mechanism  consists  of  pyrolytic 
conversion  of  solid  reactant  to  gas-phase  products.  The  other  is  that  the 


ZJ.D.  Buckmaster ,  A.K.  Kapila  and  G.S.S.  Ludfordy  "Linear  Condensate 
Deflagration  for  Large  Activation  Energy >"  Acta  Astronautica>  Vol. 
p.  593  (1976). 


14 


conversion  is  a  phase  change  between  condensed  and  gaseous  states,  i.e.,  either 
sublimation  or  evaporation  from  a  melt  layer. 

The  pyrolysis  hypothesis  is  usually  expressed  by  the  relation 


if  -E  /RT 

M  =  M  e  s  s  •  (19) 

o 

4  5 

There  have  been  few  attempts  to  estimate  M0  and  Es  in  an  a  phAjOKu  sense. 4  5 6 
The  calculation  of  Wilfong,  Penner  and  Daniels^  was  based  on  the  energy  needed 
to  break  the  0-N02  bond  in  nitrocellulose  (NC) ,  the  surface  density  of  those 
bonds,  and  application  of  transition  state  theory  to  obtain  a  frequency  factor. 
The  value  of  this  frequency  factor,  although  typical  for  gas-phase  unimolecular 
reactions,  is  five  or  more  orders  of  magnitude  lower  than  that  derived  from  NC 
thermal  decomposition  experiments. 1*^,7  m0st  cases,  however,  M0  and  Es  are 
simply  used  as  fitting  parameters  to  calibrate  a  given  burning  rate  model  to 
experimental  data. 

Evaporative  mechanisms  usually  take  the  integrated  Clausius-Clapyron  form 
for  the  propellant  vapor  pressure  (P^)  immediately  adjacent  to  the  surface, 
i .  e. , 


P 


B 


P 

o 


-E  /RT 
e  v  s 


(20) 


where  Ev  is  the  heat  of  vaporization  and  P0  is  a  constant.  Use  of  this  equili¬ 
brium-based  expression  may  have  validity  as  an  idealization  in  spite  of  the 
inherently  non-equilibrium  phenomenon  of  burning.  One  can  expect  that  equili¬ 
brium  would  be  established  at  a  rate  which  depends  on  the  molecular  velocity 
(of  molecules  escaping  the  surface).  The  vapors  outside  the  surface  are  removed 
at  a  rate  depending  on  the  transport  velocity  (V^  in  Eqn.  (2)).  The  molecular 


4R.E.  Wilfong,  S.S.  Penner  and  F.  Daniels ,  "An  Hypothesis  for  Propellant 
Burning, ,r  J .  Phys.  &  Cottoid  Chem. ,  Vot .  54,  p.  863  (1950). 


5R.D.  Schultz  and  A.O.  Dekker ,  "The  Absolute  Thermal  Decomposition  Rates  of 
Solids , "  Fifth  Symposium  (International)  on  Combustion,  The  Combustion 
Institute  (1960),  p.  260. 

6R.D.  Smith,  "Pyrolysis  of  Dissolved  Nitrocellulose,"  Rature,  Vol.  170, 
p.  844  (1952). 

?R.W.  Phillips,  C.A.  Orlick  and  R.  Steinberger,  "The  Kinetics  of  the  Thermal 
Decomposition  of  Nitrocellulose,"  J.  Phys.  Chem. ,  Vol.  59,  p.  1034  (1959). 
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velocity,  however,  should  be  fast  relative  to  the  transport  velocity  or  the 
notion  of  local  temperature  would  break  down.  It  is  reasonable,  therefore,  to 
assume  that  at  least  a  condition  of  "near"  equilibrium  vapor  pressure  might 
well  prevail  very  close  to  the  surface. 

E.  Solution  Matching  at  the  Interface 

The  conservation  laws  which  are  applied  in  obtaining  <f>G  and  <(>5,  viz., 
conservation  of  mass,  energy  and  species,  must  also  be  applied  to  the  solid- 
gas  interface.  In  addition,  the  solid  decomposition  mechanism  must  be  taken 
into  account.  The  mass  conservation  principle  is  satisfied  simply  by  requiring 
that  M  have  the  same  value  in  both  the  solid  and  gas  phases.  Energy  conser¬ 
vation  is  achieved  by  equating  the  energy  fluxes  across  the  interface,  i.e.. 


♦SCTS,M)  =  4>g CTg ,M) . 


(21) 


In  physical  terms,  this  means  that  the  energy  flux  required  by  the  solid  to 
maintain  a  steady-state  mass  regression  rate  M  at  a  surface  temperature  Ts 
must  be  supplied  by  the  heat  flux  which  the  gas-phase  processes  are  capable 
of  providing  for  the  same  values  of  M  and  Ts. 

Species  conservation  across  the  interface  is  assured  by  the  formulation 
of  an  appropriate  boundary  condition.  Since  we  are  assuming  that  no  diffusion 
occurs  in  the  condensed  phase,  the  flux  of  B  at  the  left  of  the  boundary  is 
Mmg-0.  This  value  is  reduced  on  the  gas  side  of  the  boundary  by  diffusion. 

The  continuity  condition  is  then 


M  -o  w  +0  n  ^mB 
^B  =  ^B  "  Dp  d ~ 1 


1  x=+o 


(22) 


Since  m  is  decreasing  in  the  gas  phase,  the  derivative  on  the  right  side  of 
Eqn.  (22)  will  be  negative,  so  that  mg-0  >  mg+o  as  expected.  Normally,  Eqn. 

(22)  would  be  used  as  a  boundary  condition  in  solving  Eqn.  (3)  in  order  to 
obtain  $  (Ts,M).  This  is  done  in  Section  IIIA  where  by  assuming  a  constant 

temperature  reaction  rate  an  analytic  solution  can  be  found.  It  would  also  be 
done  in  a  numerical  solution  to  Eqns.  (3)  and  (4).  In  one  of  the  approximations 
to  be  discussed  (Section  IIIC)  we  shall  find  it  convenient  to  view  Eqn.  (22) 
simply  as  a  constraint  to  be  observed  in  matching  the  heat  fluxes  (Eqn.  (21)). 

In  this  connection  it  is  useful  to  place  Eqn.  (22)  in  a  slightly  different  form. 
Using  Eqn.  (15)  to  relate  the  derivatives  of  m  and  T,  and  substituting  for 
Dp  from  Eqn.  (12),  we  find 


m 


B 


-o 


+o 

mB  +  QgM 


(23) 


In  addition  to  these  continuity  requirements  an  additional  constraint  is 
imposed  by  the  surface  decomposition  mechanism.  For  pyrolysis  M(TS)  is  given 
by  Eqn.  (19).  The  phase  change  mechanism  quantified  by  Eqn.  (20)  does  not 
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explicitly  relate  the  unknowns  TS,M,  and  m^*0,  so  that  the  nature  of  the  con¬ 
straint  is  not  clear.  Starting  from  the  relation  P  =  Pg  +  Pq,  where  Pq  is  the 
partial  pressure  of  C,  and  using  Eqn.  (10),  one  can  arrive  at  the  expression 


m.+°(T  ) 
s 


wv 


B 


p  ♦  CN2-1)PBCTS) 


(24) 


Eqn.  (24),  where  Pb(Ts)  is  given  by  Eqn.  (20),  is  then  the  most  convenient 
operational  constraint  corresponding  to  the  phase  change  mechanism. 

In  summary,  one  can  view  the  problem  of  finding  the  steady-state  solution 
as  one  in  which  the  three  unknowns,  Ts,  M,  mg+0  are  to  be  found  from  Eqns.  (21), 
(23),  and  either  Eqn.  (19)  or  Eqn.  (24).  The  way  in  which  this  is  done  will 
depend  in  detail  on  the  approximation  method  used,  a  matter  which  will  be  dis¬ 
cussed  in  Section  III. 


F .  Simplifications 

For  purposes  of  making  first  order  comparisons  of  different  theories  for 
$q,  much  of  the  detail  of  previous  sections  is  not  needed.  In  particular,  we 
wish  to  simplify  the  problem  by  holding  constant  those  variables  whose 
expected  degree  of  change  would  be  unlikely  to  seriously  influence  results. 
The  following  items  are  so  identified  and  appropriately  approximated  as 

W  =  constant 


C  = 
s 


(25) 


It  is  assumed  that  in  applying  these  approximations  one  has  sufficient  infor¬ 
mation  on  a  propellant  to  make  reasonable  estimates  for  these  constants. 

A  further  simplification  we  shall  wish  to  make  is  to  assume  negligible 
species  transport  by  diffusion.  The  quality  of  this  approximation  is  by  no 
means  self-evident,  and  can  best  be  assessed  by  comparison  to  a  calculation 
which  does  not  neglect  diffusion.  Rice  and  Ginell^  have  done  this  and  found 
that  the  character  of  the  burning  rate  behavior  is  not  seriously  altered  by 
this  assumption.  It  is  important  to  note  that  in  the  zero  diffusion  limit 
we  can  no  longer  assume  that  the  Lewis  number  is  unity  as  this  would  force 
a  neglect  of  thermal  conduction  (by  Eqn.  (12))  as  well.  With  no  diffusion 
the  species  transport  equation  (Eqn.  (3))  can  then  be  written 


dm 

M  ar ■  - R(x) 


(26) 


8O.K.  Rice  and  R.  Ginell,  "The  Theory  of  the  Burning  of  Double-Base  Rocket 
Powders , "  J.  Rhys.  &  Colloid  Chem .,  Vol.  54,  p.  885  (1950). 
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and  the  boundary  condition  of  Eqn.  (22)  becomes 


m 


+o 

B 


(27) 


For  both  pyrolysis  and  evaporation  mechanisms  this  assumption  has  the  effect 
of  reducing  the  unknowns  to  only  and  M,  eliminating  Eqn.  (23). 

For  convenience  we  label  the  approximations  expressed  in  Eqns.  (25) -(27) 
as  Level  II.  A  Level  I  calculation  is  taken  to  imply  all  of  the  assumptions 
stated  in  the  development  prior  to  Eqn.  (25).  Thus  a  Level  I  calculation  is 
our  most  general  formulation  of  the  basic  idealization.  A  Level  II  calculation 
is  an  approximate  formulation  from  which  general  trends  are  still  expected  to 
emerge.  In  this  paper  we  shall  devote  the  bulk  of  our  attention  to  the  Level 
II  approximation  since  here  our  primary  interest  is  in  qualitative  trends 
implied  by  the  chosen  idealization.  The  problem  of  applying  these  concepts 
to  actual  propellant  formulations  in  a  Level  I  calculation  requires  the  dif¬ 
ficult  choice  of  kinetic  constants  and  will  be  addressed  in  a  future  article. 


III.  GAS-PHASE  APPROXIMATIONS 

In  this  section  we  shall  consider  a  number  of  different  approaches  to 
estimating  the  gas-phase  heat  feedback  capability  <J>q.  The  utility  of  the 
form  of  Eq.  (6)  will  become  clear  since  all  of  the  approximations  discussed 
can  be  reduced  to  assumptions  as  to  the  behavior  of  q(x). 

A.  Constant  Temperature  Reaction  Rate  (CTRR) 

It  is  clear  from  Eqn.  (6)  that  if  an  adequate  approximation  for  the  heat 
release  function  q(x)  can  be  found,  (|>  can  be  calculated  for  a  given  value  of 
M.  Eqn.  (11)  shows  that  in  order  to  compute  q(x),  both  T(x)  and  mg(x)  must  be 
known.  These  two  quantities  are  interrelated  by  the  energy  conservation 
equation  (Eqn.  (4))  and  the  species  continuity  equation  (Eqn.  (3)).  If  we 
assume  that  the  reaction  proceeds  as  though  it  were  at  some  constant  temperature 
T  (historically  taken  as  the  flame  temperature,  T^)  and  with  constant  W,  i.e.. 


R(x)  =  Ro  m^ 


where 


(28) 


R  5 


WvPVAn 


-eg/rti 


RVT^ 


(29) 


then  Eqn.  (3)  can  be  solved  exactly  for  v=l  (using  Eqn.  (22)  as  interface 
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boundary  condition) . 


The  solution  is 


m 


-o  -ax 
mB  6 


B  ,,  Dp  . 


(for  v=l) 


(30) 


where 


Using  Eqns.  (28)  and  (30)  we  may  then  arrive  at  an  expression  for  the  gas- 
phase  heat  feedback  by  integrating  Eqn.  (6). 


<tv  = 


X(^GRomB 


Dpa.  Xa 
MC  (1+  -—•)  (1+  Tjr^- 
p  M  MC 


(for  v=l) 


(31) 


It  is  of  interest  to  compare  this  solution  which  includes  the  effects  of 
species  diffusion  to  one  in  which  such  effects  are  neglected.  This  can  be 
accomplished  by  solving  Eqn.  (3)  without  the  second  derivative  term  (or  taking 
the  limit  of  Eqn.  (30)  as  Dp  0).  Again  for  v=l,  the  results  are 


m 


B 


(v=l,  no  diffusion) 


(32) 


and 


X(^GRomB 


XR 

MCp(1+  -T2) 
P  MC 

P 


(v=l,  no  diffusion). 


(33) 


The  functional  forms  of  Eqns.  (31)  and  (33)  bear  some  resemblance,  and  in  the 
limit  of  small  R0/M 2  the  forms  become  identical.  Physically,  this  limit  corre¬ 
sponds  to  the  heat  feedback  being  small  compared  to  the  maximum  possible 
heat  feedback  Q^Mnig-0.  (See  Section  IVA1)  . 

For  the  case  of  a  second  order  reaction,  i.e.,  v=2,  Eqn.  (3)  appears  not 
to  have  a  closed  form  solution,  However,  it  can  be  solved  analytically  in  the 


19 


limit  of  zero  diffusion,  in  which  case 


-o 

(v=2,  no  diffusion) #  (34) 


for  <|>  is 


E2(^-CP  _^)  (v=2,  no  diffusion)  (35) 

o  B 

exponential  integral. 

Equation  (35)  can  be  related  to  two  of  the  early  attempts  to  model  homo¬ 
geneous  solid  propellant  combustion.  Both  Parr  and  Crawford9  and  Rice  and 
Ginell8  assumed  a  second  order  reaction  in  the  "fizz  zone"  or  gaseous  region 
immediately  adjacent  to  the  surface  of  a  burning  double-base  propellant.  Parr 
and  Crawford9  assumed  zero  diffusion  and  arrived  at  a  first  approximation 
(Eqn .  (42)  of  Ref.  9)  which  is  essentially  the  same  expression  as  a  mixture  of 
our  Eqns .  (18)  and  (35).  Rice  and  Ginell8  at  first  treated  the  zero  diffusion 
case  and  also  arrived  at  an  equivalent  of  our  Eqns.  (18)  and  (35)  (Cf.  Eqns. 

(9),  (10)  and  (11a)  of  Ref.  8).  These  latter  authors  then  introduced  a 
diffusion  correction  for  the  second  order  reaction  case  and  found,  as  have  we 
for  v=l  above,  that  the  functional  forms  of  the  solution  with  and  without 
diffusion  are  similar.  Rice  and  Ginell8  justified  the  use  of  a  constant 
temperature  reaction  rate  as  a  limiting  case  of  zero  activation  energy  (Eq) 
and  simply  used  kG  as  an  adjustable  constant  in  fitting  their  theory  to  burning 
rate  data.  Parr  and  Crawford,9  while  conceding  that  zero  activation  energy  is 
the  proper  limiting  case  corresponding  to  use  of  Eqn.  (29),  chose  to  retain 
explicitly  both  the  activation  energy  (Eq)  and  frequency  factor  (Aq)  as  adjust¬ 
able  parameters  in  fitting  their  theory  to  burning  rate  data.  Their  rationale 
was  to  consider  the  constant  temperature  reaction  rate  as  an  approximate  for¬ 
mulation.  In  both  treatments8*9  the  temperature  Ti  was  taken  to  be  that  at 
the  end  of  the  fizz  zone,  i.e.,  in  our  notation. 

B.  Delta  Function  Heat  Release  (DFHR) 


D 

mB  =  m  -°R  x 
B  o 

O  - M-> 


The  corresponding  expression 

M2Cp 
XR  m-° 

_  -o  o  B 

*G  =  V"B  e 


where  (z)  is  the  standard 


Another  approach  to  evaluating  <|>q  is  to  assume  that  the  heat  release 
function  q(x)  can  be  approximated  by  a  Dirac  delta  function,  i.e., 


qO)  =  qQ  <s(x-x*) 


(36) 


9R.G.  Parr  and  B.L.  Crawford ,  Jr.,  "A  Physical  Theory  of  Burning  of  Double- 
Base  Pocket  Rocket  Propellants,"  J.  Phys.  &  Colloid  Chem. ,  Vol.  54,  p.  929, 
1950. 
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where  x*  is  the  so  called  flame-standoff  distance  and  q  can  be  evaluated  by 
.  o 

the  requirement 


f  q(x)dx  =  QgMmB 
J  o 


(37) 


which  is  simply  the  flux  corresponding  to  the  total  gas-phase  heat  release. 
Eqn.  (37)  can  be  obtained  by  integrating  Eqn.  (3)  over  (+0,°°)  subject  to  the 
boundary  condition  Eqn.  (22).  It  follows  that  the  constant  in  Eqn.  (36)  is 
given  by 


%  =  ‘ 


(38) 


Equation  (36)  is  a  reasonable  picture  of  q(x)  in  the  limit  of  infinitely  large 
activation  energy  E q,  as  can  be  seen  graphically  in  the  laminar  flame  theory 
paper  of  Spalding. 1™  The  difficulty  of  using  Eqn.  (36)  is  not  in  its  functional 
form  but  in  evaluating  a  proper  value  for  x*.  Inserting  Eqn.  (36)  into  Eqn. 

(6)  we  obtain 


MC 

_ E  x* 

4>g  =  Qg^b"0  e  ^  *  (39) 


Thus  the  value  of  x*  is  critical  to  this  approximation. 

As  representative  of  the  "f lame-sheet,T  models  we  shall  discuss  only  the 
model  devised  by  Hermancell  and  later  incorporated  into  the  BDP  composite 
propellant  model. 12 >13  All  of  these  models  share  the  characteristic  that  the 
reaction  rate  is  evaluated 
in  detail  sufficiently  that 

BDpl^l^  formulation.  The  other  flame-sheet  models,  of  which  this  writer  is 


10 

D.B.  Spalding,  "II.  One-Dimensional  Laminar  Flame  Theory  for  Temperature- 
Explicit  Reaction  Rates,"  Combustion  &  Flame,  Vol.  1,  p.  296  (1957). 

11 

C.E.  Hermance,  "A  Model  of  Composite  Propellant  Combustion  Including 
Surface  Heterogeneity  and  Heat  Generation,"  AIAA  Journal,  Vol.  4,  p.  1629 
(1966). 


only  at  the  flame  temperature;  however,  they  differ 
our  results  will  pertain  only  to  the  Hermancell/ 


12 

M.W .  Beckstead,  R.L.  Derr  and  C.E .  Price,  "A  Model  of  Composite  Solid- 
Propellant  Conibustion  Based  on  Multiple  Flames,  "  AIAA  Journal,  Vol.  8, 
p.  220  (1970). 


13 

M.W.  Beckstead ,  R.L.  Derr  and  C.E.  Price 3  "The  Combustion  of  Solid  Mono¬ 
propellants  and  Composite  Propellants , "  Thirteenth  Symposium  (International) 
on  Combustion ,  The  Combustion  Institute  (1971),  p.  1047. 
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aware,  are  due  to  Culick^  and  Williams.  In  Culick's  model,  q  is  not 
determined  by  Eqn.  (37)  but  specified  a  pfUofii  as  QqRq  (defined  in°Eqn.  29  with 
Tl  =  Tf) .  The  argument  of  the  delta  function  is  then  expressed  in  a  trans¬ 
formed  coordinate.  With  this  specification  of  q0,  however,  one  can  show  that 
the  derived  burning  rate  formula  depends  on  the  coordinate  transformation 
used,  thereby  invalidating  the  method.  The  Williams  model^  avoids  the 
delta  function  formalism  by  applying  asymptotic  analysis  to  the  flame  sheet 
idealization.  The  resulting  solution  is  functionally  similar  (though  not 
identical)  to  the  Hermance  formulation. 

Hermance^  argues  for  the  value  of  x*  as  follows. 


x*  =  Ut 


(40) 


where  U  is  the  average  velocity  of  gases  flowing  away  from  the  solid  interface 
and  t  is  the  average  reaction  delay  time.  These  quantities  are  in  turn  esti¬ 
mated  as 


_  M_  (41) 

"  pf 


where  is  the  gas  density  at  the  flame  temperature,  and 


T 


1 


e 


-EG/RTf  ' 


(42) 


Combining  (40),  (41)  and  (42)  we  obtain 


M 


v.  -E„/RT_ 
p  e  G  f 


f  G 


(43) 


Hermance  considered  only  v=2;  therefore  we  have  generalized  his  result  to 
arbitrary  reaction  order.  Use  of  Eqn.  (42)  would  seem  to  be  inconsistent 
with  the  sense  in  which  the  delta  function  might  be  considered  a  valid 
idealization,  i.e.,  in  the  limit  of  high  Eq.  In  this  limit  the  reaction  is 
essentially  confined  to  a  thin  spatial  zone  whose  thickness  is  small  compared 
with  its  distance  from  the  surface.  The  t  required  by  Eqn.  (40)  is  the  mean 

14F.E.C.  Culick ,  "An  Elementary  Calculation  of  the  Combustion  of  Solid 
Propellants,"  Astronautica  Acta,  Vol.  14,  p.  171  (1969). 

1 5F.A .  Williams,  "Quasi- Steady  Gas-Phase  Flame  Theory  in  Unsteady  Burning 
of  a  Homogeneous  Solid  Propellant,"  AIAA  Journal,  Vol.  11,  p.  1S28  (1973). 
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time  it  takes  for  the  gaseous  reactants  to  travel  from  the  surface  to  the 
thin  reaction  zone.  The  t  represented  by  Eqn.  (42],  on  the  other  hand,  is 
related  to  the  time  it  takes  to  complete  the  reaction  within  the  reaction 
zone.  This  same  expression  is  used  by  Rice  and  Ginell8  to  describe  the  zero 
activation  energy  limit.  It  would  seem,  therefore,  that  the  true  value  of 
x*  is  underestimated  by  Eqn.  (43],  leading  to  an  overestimate  of  <(>£.  In  the 
more  rigorous  treatment  of  the  high  EG  limit  by  Williams15,  however,  an 
expression  for  x*  can  be  obtained  which  is  the  same  as  Eqn.  (43)  except  for 
a  multiplicative  constant  (more  precisely,  a  quantity  which  is  nearly  constant) 
Such  a  circumstance  seems  to  justify  use  of  Eqn.  (43);  however,  in  an  actual 
numerical  calculation  the  multiplicative  constant  may  have  a  significant 
effect  since  x*  appears  in  an  exponent.  As  we  are  primarily  interested  in 
functional  relationships  in  this  paper  we  may  ignore  such  fine  points  and 
treat  the  Hermance  flame  standoff  as  adequately  representative  of  the  delta 
function  approach.  Substituting  Eqn.  (43)  into  Eqn.  (39),  one  obtains 


$G  =  Qg^B  °  exP 


M2Cp  RVT^  ^ 
XWVPVAg  e  "EG//RTf  J 


(44) 


C.  Quasi-Constant  Heat  Release  (QCHR) 


We  have  developed  an  approximation  to  which  is  complementary  to  the 
CTRR  and  DFHR  limiting  cases.  The  spatial  variation  of  the  heat  release 
function  is  essentially  determined  by  a  competition  between  the  decreasing 
reactant  mass  fraction  mg  and  the  increasing  reaction  rate  kQ(cf.  Eqn.  (7)). 

In  the  CTRR  case  the  mg  variation  dominates  so  the  q(x)  is  monotonic  decreasing 
In  the  DFHR  model  kQ  initially  increases  faster  than  m^  decreases  so  that 
q(x)  increases  over  most  of  the  distance  in  the  interval  (0,x*).  Ultimately, 
however,  Tf  is  approached  and  mg  drives  q(x)  to  zero.  Our  approximation 
assumes  that  the  rate  of  increase  of  kg  is  matched  by  the  rate  of  decrease 
of  nig  at  least  out  to  some  distance  x^  >>  ^/MCp,  i.e.,  q(x)  is  assumed  constant 
from  the  surface  to  X]^.  Equation  (6)  can  then  be  written 


q(x)  dx 


i  Xq(o) 
MC 

P 


e 


xq(o) 

MC 

P 


q(x1+y)  dy 


(45) 


where  y  =  x  -  Xj.  Since  we  have  assumed  that  x^>>  X/MCp,  the  second  term  in 
Eqn.  (45)  will  be  small  due  to  the  exponential  and  can  be  neglected  (assuming 
also  that  q(x)  does  not  increase  too  rapidly  for  x  >  x^) .  Eqn.  (45)  then 
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becomes 


*G 


*q(o) 

MC 


p 


(46) 


This  formalism  is  the  same  as  that  used  in  the  steady-state  component  of 
the  KTSS  non-steady  combustion  theory. ^  However,  in  that  theory  the  equiva¬ 
lent  of  Eqn.  (46)  was  not  used  to  calculate  a  pnjjOhA,  the  steady-state  burning 
behavior.  Instead,  the  pressure  dependence  of  the  term  Xq(o)/Cp  was  deduced 
semi-empirically  by  imposition  of  an  assumed  burning  rate  law  of  the  form  M=apn. 

Equation  (46)  also  bears  a  formal  resemblence  to  a  simple  burning  rate 
theory  developed  by  Kubota,  et  al.-^  as  a  baseline  in  their  study  of  double 
base  catalysis.  In  that  work,  however,  q(x)  was  assumed  constant  throughout 
the  flame  zone.  The  constant  value  was  then  evaluated,  not  at  the  interface, 
but  at  the  flame  temperature  and  density.  In  order  to  avoid  the  result  that 
q(Tf)  =  0,  because  the  reactant  mass  fraction  is  zero  at  the  flame  temperature, 
the  mass  fraction  appears  to  have  been  arbitrarily  reset  to  its  initial  value. 
The  net  result  is  a  model  which  is  identical  to  the  low  heat-feedback  limit 
of  the  (zero  diffusion)  CTRR  formulation  (see  Section  IVA1). 

Evaluating  Eqn.  (11)  at  x  =  +o  and  substituting  the  result  into  Eqn.  (46) 
gives  <J>g  as  an  explicit  function  of  Ts,  M,  and  mg+^.  This  enables  a  Level  I 
solution  to  be  found  by  simultaneously  solving  Eqns.  (21),  (23)  and  either 
Eqn.  (19)  or  (24),  depending  on  the  desired  surface  decomposition  mechanism, 
for  the  three  unknowns  Ts,  M,  and  mB+0.  Thus,  the  net  effect  of  the  QCHR 
approximation  is  to  replace  a  difficult  non-linear  differential  equation  problem 
with  a  relatively  simple  algebraic  one. 

A  Level  II  formulation  (see  Section  IIF)  of  the  QCHR  approximation  can 
be  shown  to  give 


-E  /RT 

XQg  (W  P  mB"°)  Ag  e 

MC  RVTV 
P  s 


(Level  II)  (47) 


Unlike  the  CTRR  and  DFHR  models  the  criterion  for  validity  of  the  QCHR 
assumption  is  not  clearly  expressible  in  terms  of  limiting  values  of  the 
input  parameters.  Rather,  the  specification  is  placed  on  the  behavior  of 
an  integrated  quantity,  q(x),  which  to  some  degree  involves  all  of  the  input 
parameters.  Since  the  kinetics  parameters  pertinent  to  burning  propellants 

16H.  Krier,  J.S.  Tien,  W.A.  Sirignano  and  M.  Summer  fie  Id,  "Nonsteady  Burning 
Phenomena  of  Solid  Propellants:  Theory  and  Experiments ,"  AIM  Journal , 

Vol.  6,  p.  278  (1968). 

17  N.  Kubota ,  T.J.  Ohlemiller ,  L.H.  Caveny  and  M.  Summerfield,  "The  Mechanism 
of  Super-Rate  Burning  of  Catalyzed  Double  Base  Propellants,"  AMS  Report 
No.  1087,  Princeton  University,  March  1973. 
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are  not  reliably  known,  this  method  of  specification  may  not  represent  a 
disadvantage.  In  fact,  observations  of  the  burning  process  inevitably  sense 
integrated  quantities,  so  that  it  may  be  easier  to  judge  the  appropriateness 
of  the  QCHR  approximation  than  those  which  specify  a  limiting  value  of  E^. 
This  approach  is  illustrated  in  the  next  section. 

D.  Review  of  Experimental  Evidence  for  Gas  Phase 


The  most  detailed  experimental  investigations  of  the  thermal  wave  struc¬ 
ture  in  solid  propellant  combustion  have  been  measurements  of  temperature 
profiles  obtained  by  imbedding  small  thermocouples  in  propellant  samples. 
Interpretration  of  these  measurements  is  beset  by  problems  such  as  thermal 
inertia,  thermal  conduction  through  the  junction  leads,  and  determining  when 
the  thermocouple  emerges  from  the  solid.  Because  of  the  steep  temperature 
gradients  actually  measured,  these  problems,  which  are  still  not  fully 
resolved,  can  present  serious  obstacles  to  extracting  accurate  temperature 
profiles  from  the  raw  data.  Nonetheless,  these  measurements  can  provide 
valuable  qualitative  and  some  quantitative  insight  as  to  the  nature  of  the 
solid  propellant  burning  process.  They  have,  for  example,  revealed  the 
existence  of  distinct  plateaus  in  the  spatial  temperature  profile  in  the 
gas  phase  for  single  and  double  base  propellants. 17, 18, 19,20  The  generally 
accepted  picture  (for  these  propellants)  has  emerged  of  two  exothermic  kinetic 
stages  —  one  immediately  adjacent  to  the  surface  (fizz  zone)  and  one 
associated  with  the  visible  flame.  These  two  zones  are  separated  by  a  "dark” 
zone  characterized  by  temperature  plateau  where  presumably  the  flame  zone 
reactants  are  created  in  thermoneutral  reactions. 

Evidence  as  to  whether  or  not  the  visible  flame  reactions  affect  the 
burning  rate  is  mixed.  One  test  is  to  look  for  a  kink  in  the  log  r  vs.  log  P 
curves  through  the  pressure  range  where  the  visible  flame  first  appears. 
Crawford,  et  al.^l  and  Heath  and  Hirst"*-  both  made  particular  note  of  the 
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absence  of  such  a  kink.  Suh  and  Clary  ,  on  the  other  hand,  found  a 
prominent  slope  change  at  the  onset  of  the  visible  flame.  All  of  the 
measurements  were  made  with  double-base  propellants  but  compositions  varied. 

21 

While  the  energy  release  in  the  flame  zone  is  substantial  ,  we  shall 
assume  that  it  contributes  negligibly  to  <J>g.  Further  we  assume  that  the  fizz 
zone  reaction  or  reactions  can  be  adequately  represented  by  our  single  gas 
phase  reaction  B+C.  Thus  is  to  be  taken  as  the  dark  zone  temperature  and 
not  the  final  propellant  flame  temperature. 

Klein,  et  al.,15  first  used  the  temperature  profiles  generated  by  imbed¬ 
ded  thermocouples  to  estimate  the  heat  release  function  q(x)  from  the  energy 
conservation  equation  (Eqn.  (4)).  His  results  for  nitrocellulose  strands  at 
28  atm  pressure  are  given  in  Fig.  2  which  shows  the  inferred  fizz  zone  energy 
release.  We  have  taken  the  thermocouple  data  of  Kubota,  et  al.,1^  for  a 
laboratory  double-base  propellant  PNC/TMETN  (particulate  nitrocellulose  and 
trimethylolethane  trinitrate)  which  has  similar  burn  rate  characteristics  as 
NC/NG  (nitrocellulose  and  nitroglycerin)  and  constructed  the  function  q(x) 
in  the  manner  of  Klein,  et  al.1®  This  function  is  computed  at  1.2  atm  and 
21  atm  and  shown  in  Figs.  3  and  4,  respectively. 

Since  the  QCHR  approximation  is  couched  in  terms  of  q(x),  the  "empirical" 
q(x)  functions  can  be  used  to  test  the  model  directly.  The  CTRR  and  DFHR 
Models  have  also  been  described  in  terms  of  q(x)  for  this  same  purpose. 

Using  Eqns .  (5),  (28),  and  (32)  or  (34)  one  can  see  that  q(x)  in  the 
CTRR  Model  is  a  monatonic  decreasing  function  of  x.  This  behavior  is  not 
well  approximated  in  Figures  2-4.  The  rising  nature  of  q(x)  results  from 
a  reaction  rate  which  is  increasing  (due  to  temperature  rise)  at  a  faster 
rate  than  the  reactant  is  being  consumed. 

The  self-consistency  of  the  DFHR  Model  can  be  partially  tested  by 
integrating  the  q(x)  functions  over  x  to  obtain  an  estimate  of  Q^Mnig0 
(Eqn.  (37)),  then  using  the  position  of  the  q(x)  peak  as  x*.  A  "model"  value 
for  <J)q  can  then  be  computed  from  Eqn.  (.39) .  This  can  be  compared  with  the 
"empirical"  value  of  <J>g  obtained  by  performing  the  integration  in  Eqn.  (6)  . 

The  results  of  such  a  calculation  are  shown  in  Table  1.  (The  error  associated 
with  these  integrations  is  on  the  order  of  10%  as  judged  by  the  self-consistency 
between  the  measured  temperature  gradient  at  the  surface  times  \  and  the 
integration  over  q(x)).  Although  the  q(x)  curves  in  Figures  2-4  are  not 
characterized  by  a  small  peak-width-to-standoff-distance  ratio,  the  "model" 
and  "empirical"  values  agree  quite  well.  Not  probed  by  this  test,  however, 
is  the  accuracy  with  which  the  actual  flame  standoff  distance  would  be 
predicted  by  the  idealized  value  as  given  by  the  Williams  model15,  for  instance. 

The  constant  q  assumption  of  the  QCHR  model  also  appears  to  be  at  odds 
with  the  figures;  however,  the  value  of  the  model  as  an  approximation  depends 
on  the  extent  to  which  q(x)  changes  over  the  distance  scale  X/MCp.  The  QCHR 
model  value  of  is  obtained  from  Figures  2-4  using  Eqn.  (46).  As  shown  in 


2SN.P.  Suh  and  D.L.  Clary ,  "Steady-State  Burning  of  Double-Base  Propellants 
at  Low  Pressures,"  AIAA  Journal,  Vol .  8,  p.  825  (1970). 
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Figure  2.  Empirical  gas-phase  heat  release  function  calculated  by 
Klein,  et  al.,'^  for  NC  at  28  atm. 
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Figure  3.  Empirical  gas-phase  heat  release  function  based  on  data  of 
Kubota,  et  al . ,  for  PNC/TMETN  at  1.2  atm. 
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Figure  4.  Empirical  gas-phase  heat  release  function  based  on  data  of 
Kubota,  et  al.J?  for  PNC/TMETN  at  21  atm. 
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TABLE  1.  COMPARISON  OF  HEAT  RELEASE  MODELS 
USING  EMPIRICAL  HEAT  RELEASE  FUNCTIONS 


u 

o 

u 

u 

PJ 

o\° 


cc 

<y 


o\° 

o\° 

©\° 

O 

O 

O 

CM 

sO 

s — / 

v—/ 

S— * 

sO 

oJ 

to 

SO 

CM 


O 

o 


CO 

PJ 

ID 

rJ 

< 

> 

us 

< 


CO 

Q 

PJ 

m 

pl 


as 


* 

o 

Vi 

Vi 

PJ 


£ 

o 

•H 

+J 

o 

£ 

£ 

Uh 

« 

rH 

o 

Q 


Oh 

e 

pj 


o\° 

to 

r>. 


o\° 


CM 

o\ 

• 

• 

to 

r^. 

oo 


o\° 

o 


to 


to 


to 

ti 

. 

*H 

oo 

CL, 

•rH 

s _ s 

CL, 

r — \ 

_ / 

o 

CM 

s 

u 

S 

£ 

• 

£ 

+-» 

oo 

£ 

in 

•H 

CM 

o 

Cl. 

• 

rH 

u 

V _ / 

rH 

CM 

cm 

6 

r 

r 

r 

+j 

r-. 

+J 

£ 

rH 

rH 

£ 

2 

2 

£ 

00 

H 

H 

i— i 

CM 

PJ 

PJ 

rH 

<D 

•s 

p 

p 

Oh 

00 

O 

rH 

U 

u 

Vl 

u 

.  Z 

2 

Cm 

2 

CL 

CL 

30 


Table  1,  the  QCHR  values  of  $q  seriously  underestimate  the  empirical  value 
in  the  low  pressure  case  (Figure  3)  but  are  quite  credible  at  the  higher 
pressures.  The  Kubota,  et  al.l^  data  did  not  extend  to  higher  pressures  but, 
taken  together,  Figures  3  and  4  might  well  represent  a  trend  that  the  QCHR 
approximation  improves  with  pressure.  Since  the  basis  of  the  approximation 
is  that  q(x)  remains  relatively  constant  for  a  distance  away  from  the  surface 
on  the  order  of  A/MCp,  it  is  equivalent  to  there  being  little  curvature  in  the 
temperature  profile  out  to  this  distance.  Thus,  if  the  temperature  rises 
linearly  away  from  the  surface  for  a  distance  on  the  order  of  A/MCp,  the  QCHR 
model  should  give  a  reasonable  approximation  to  <|>  .  The  thermocouple  data  of 
Zenin^  for  a  NC/NG  propellant  (ballistite)  give  similar  results  as  that  of 
Kubota,  et  al.,17  at  20  atm  pressure.  Zenin^O  also  obtained  data  at  pressures 
up  to  150  atm.  At  50  atm  the  temperature  rise  appears  to  be  quite  linear  to 
a  distance  well  beyond  A/MCp,  tending  to  support  the  prospect  of  QCHR  improving 
with  pressure.  At  higher  pressures  the  resolution  is  too  coarse  to  judge  the 
linearity  of  T. 

In  summary,  the  available  thermocouple  data  gives  little  encouragement  to 
the  CTRR  approximation  since  the  shape  of  q(x)  over  the  region  of  influence  is 
incorrectly  predicted.  The  delta  function  formulation  appears  to  be  a  reason¬ 
able  idealization  of  q(x),  but  we  were  unable  to  test  the  adequacy  of  the 
idealized  flame  standoff  distance. The  QCHR  approximation  appears  to  be 
reasonably  consistent  with  the  thermocouple  data  above  ^  20  atm  for  nitrate 
ester  propellants.  We  caution,  however,  that  comparisons  with  this  data  are 
made  in  part  to  demonstrate  what  we  feel  to  be  a  promising  analysis  procedure. 
We  have  made  no  attempt  to  assess  the  reliability  of  the  embedded  thermocouple 
method  for  this  purpose. 


IV.  ANALYTIC  BURNING  RATE  EXPRESSIONS 

In  this  section  we  shall  mate  the  various  approximations  for  <|>g  with  $$ 
for  different  surface  decomposition  models  in  order  to  derive  expressions  for 
the  burning  rate.  We  are  interested  foremost  in  trying  to  determine  the  theo¬ 
retical  basis  of  the  ubiquitous  empirical  burning  rate  laws  which  have  been 
variously  determined^ >  21  as 


r  =  a  +  bP 


(48) 


r  =  bPn  (49) 

r  =  a  +  bPn  (50) 


where  a,  b,  and  n  are  constants  determined  from  fits  of  these  functional  forms 
to  the  burning  rate.  In  addition  we  shall  be  interested  in  the  temperature 
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sensitivity  of  the  burning  rate  at  constant  pressure  which  is  defined  by 


a 

P 


P 


(51) 


where  the  subscript  p  denotes  that  the  derivative  is  taken  while  holding  the 
pressure  constant.  The  most  common  empirical  representation  of  a  is1^4 


a 

P 


(52) 


where  the  constants  c  and  d  are  fitting  parameters.  In  order  to  maintain  a 
maximum  of  analytic  flexibility,  we  shall  confine  ourselves  to  estimates  based 
on  Level  II  approximations  in  every  case. 

A.  Pyrolysis  Decomposition  Mechanism 

As  discussed  in  Section  IIE  one  may  obtain  a  value  for  the  burning  rate 

at  some  (T0,  P)  by  equating  <j>g  (Eqn.  18)  to  the  model  of  <|>g  of  interest  (Eqns.  33, 

35,  44,  or  47).  Using  the  pyrolysis  law  (Eqn.  19)  to  eliminate  M  as  an  independ¬ 
ent  variable,  <|>s  can  be  expressed  as  a  function  of  Ts  and  T0  only.  Any  of  the 

models  for  can  be  expressed  as  a  function  of  Ts  and  P  only.  For  the  CTRR  and 
DFHR  models  in  which  depends  explicitly  on  Tf,  Tf  may  be  found  as  a  function 
of  Ts  by  solving  the  equation 

Tf  -  Ts  *  <Vb'0  -  WTf>/M)/Ep 


which  equates  the  heat  retained  in  the  gas  to  the  difference  between  the  total 
gas-phase  heat  release  and  the  heat  lost  from  the  gas  to  the  surface  by 
conduction.  One  then  must  find  that  value  of  Ts  for  which  <f>s(Ts,T0)  =  CTS  P) 
and  use  the  pyrolysis  law  to  obtain  the  burning  rate. 

It  is  informative  to  illustrate  this  matching  process  graphically  by 
plotting  both  <J>g  and  <J>g  as  a  function  of  Ts.  This  is  done  in  a  schematic 
fashion  in  Figure  5,  where  the  intersection  of  the  <{>  and  <f>Q  curves  occurs  at 
the  surface  temperature  pertinent  to  either  (T0,  Pi)  or  (T0,  P2)  •  Note  that 
this  construct  effects  a  separation  of  the  influences  of  pressure  and  initial 
temperature  on  the  problem  since  <(>„  does  not  depend  on  T0  and  <f>s  does  not 
depend  on  P.  For  values  of  Ts  not  close  to  (where  =  T0  +  Qs/Cp  is  the 
surface  temperature  for  which  no  gas-phase  heat  feedback  is  required! ,  4>S  i-s 
dominated  by  the  exponential  character  of  M  (by  Eqn.  (19));  but  for  Ts 


24K.K.  Andreyev  and  A.F.  Belyayev,  Theory  of  Explosive  Substanoesj  Moscow 
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log  [+,«T,.y] 


'$ 

SURFACE  TEMPERATURE 


Figure  5.  Illustration  of  matching  condensed  and  gaseous  phase  heat  fluxes 
in  order  to  obtain  surface  temperature  for  a  given  (T0,P)  in 
the  case  of  a  pyrolysis  mechanism. 
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the  (Ts  -  T|)  terra  dominates,  driving  <j>s  to  zero  as  Ts  -*■  T^.  The  dependence 
of  <j>G  on  P  is  less  complicated  in  that  varying  P  serves  only  to  translate 
the  <j>G(Ts)  curve  vertically  with  respect  to  <|>S(TS).  To  illustrate  this,  two 
<j)G  curves  are  drawn  in  Figure  5  for  the  different  pressures  P^  and  P2 
(where  P^  >  P2) -  In  the  limit  of  sufficiently  low  pressures,  the  intersection 
of  <j>G  and  <j>s  occurs  at  an  essentially  constant  value  of  Ts  =  T^.  Thus  in  the 
low  pressure  limit  M  -*■  constant,  that  constant  depending  only  upon  M0  and  T^ 
and  not  on  the  gas-phase  parameters.  In  this  limit,  it  is  fair  to  say  that  the 
solid  reaction  is  rate  limiting.  However,  at  high  pressures  (where  M  dominates 
<f>S)  neither  the  gas  nor  the  solid  phase  reaction  alone  can  be  said  to  be  rate 
limiting,  i.e.,  M  will  depend  on  changes  in  either  the  solid  or  gas-phase 
reaction  rates  (or  heats) . 

1.  CTRR  Approximation  to  <j> 

G 

In  order  to  discuss  the  burning  rate  derived  from  the  CTRR  approximation 
to  4>  ,  we  first  show  that  Eqns.  (33)  and  (35)  can  be  placed  in  a  similar 
form.  We  introduce  the  following  analytic  approximation  for  the  exponential 
integral 

-z 
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1  +  z 


E2(z)  * 


This  approximation  gives  the  correct  limiting  values  for  E2  at  z=0  and  z-*» 
and  overestimates  E2  by  a  maximum  of  25%  at  x=0.5.  Using  this  expression  to 
reduce  Eqn.  (35),  one  can  show  that  for  v=l  or  2, 
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where  R0  is  defined  in  Eqn.  (29).  Equating  this  <j>  to  <j>  from  Eqn.  (18), 
solving  for  M  produces  the  expression 
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For  most  propellants  it  can  be  argued  that  ES/RTS  >>1,  in  which  case  the  , 
variation  of  exp  (-ES/RTS)  with  Ts  is  enormous  compared  with  that  of  (Ts-Ts  ) 
in  the  high  pressure  limit  where  TS^TS'.  The  pressure  dependence  of  the  burn 
rate  M  is  therefore  seen  to  be  approximately  M  <*  pv/^,  a  dependence  often 
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associated  with  simple  propellant  combustion  models.  ’  ^ ^ ^  This  pVy^ 
dependence  is  a  fairly  rigorous  result27  for  a  premixed  laminar  flame  with  a 
reaction  of  order  v,  though  the  boundary  conditions  in  that  case  are 
different  than  for  a  solid  propellant.  Combining  the  low  and  high  pressure 
limiting  behaviors,  M  as  calculated  in  the  CTRR  approximation  for  a  pyrolytic 
mechanism  has  approximately  the  same  overall  pressure  dependence  as  the  form 


M  = 


+  b  P 
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and  b^  decreases  slowly  with  increasing  pressure. 

The  temperature  sensitivity  in  the  high  pressure  limit  can  be  obtained  as 
follows . 
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Eqn.  (19)  may  be  used  to  obtain  the  derivative  of  M  with  respect  to  T?. 
Eliminating  M  between  Eqns.  (19)  and  (54)  gives  an  expression  containing  the 
variables  Ts,  T0,  and  Tf.  Tf  can  be  written  as  a  function  of  T0  in  Level  I 
generality  by  solving  the  equation  h^(T0)  =  h^fTf).  For  the  case  of  constant 
and  equal  specific  heats,  the  result  is 


Qt 

Tf  =  Tq  +  (55) 

CP 


where  Qy  =  (Qs  +  mg°  Qq)  is  the  total  heat  released  in  both  the  solid  and  gas 
phases  (through  the  fizz  zone,  for  example)  per  unit  mass  of  propellant. 

Using  this  relation  to  eliminate  Tf  as  well,  one  obtains  an  equation  involving 
only  Ts  and  T0.  Taking  the  natural  logarithm  of  both  sides  of  this  equation 
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then  performing  an  implicit  differentiation  with  respect  to  Tq  yields 
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The  parameter  3  is  related  to  the  fractional  heat  feedback  from  the  gas  phase, 
being  large  when  this  fraction  is  large  and  close  to  unity  in  the  limit  of  low 
fractional  heat  feedback.  Two  limits  can  be  identified.  For  high  fractional 
heat  feedback,  i.e.. 


*G 


*  1, 


^  1 


and  a  in  this  limit  can  be  shown  to  be 
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(56) 


a 
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We  do  not  expect  this  limit  to  be  significant  in  a  practical  sense,  however, 
since  if  <f)(j  =  Q^Mmg  ,  Ts  would  approach  a  constant  (as  can  be  seen  by  using 
Eqn.  (18))  and  M  would  therefore  lose  its  pressure  dependence.  Experimentally, 
most  homogeneous  propellants  show  no  sign  of  departing  from  M  «  pn  at  high 
pressure.  For  low  to  moderate  fractional  heat  feedback,  and  ES/RTS  >>1, 
the  leading  terms  in  Op  become 
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The  first  term  on  the  right  hand  side  of  Eqn.  (57)  is  similar  to  the  empirical 
form  of  Eqn.  (52).  However,  while  the  third  term  on  the  RHS  of  Eqn.  (57)  can 
likely  be  neglected  compared  with  the  second  term,  the  second  term  cannot  be 
expected  to  be  negligible, 
to  be  -V  1500°K17>T8’20  and 

reactions^, 33)  ?  Eq/2RT£2  = 

also  that  the  second  term  decreases  with  increasing  T0  contrary  to  the^empirical 
behavior  in  Eqn.  (52). 

2.  DFHR  Approximation  to 

Equating  Eqns.  (18)  and  (44)  and  solving  for  M  produces  the  following 
result  applicable  to  the  high  pressure  limit  (i.e.,  where  TS^TST).  The  term 
in  brackets  can  be  expected  to  decrease  slowly  with  pressure. 


For  example,  for  a  NC-based  propellant,  taking  Tf 
Eq  as  15  kcal/mole  (representative  of  N02-aldenyde 
.0017,  which  is  of  the  same  magnitude  as  oT 
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Again  at  low  pressures  M  approaches  the  same  constant  (a^)  as  was  found  in 
Section  IVA1. 


For  the  6 -function  case  one  finds  that 
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Like  g  of  the  preceding  section,  g  is  related  to  the  fractional  heat  feedback 
from  the  gas  phase  and  is  of  order  unity  as  long  as  is  not  very  near  Q^Mm^-0. 
We  previously  argued  that  this  circumstance  is  unlikely  to  be  encountered  over 
the  conditions  of  practical  interest.  Therefore  in  the  high  pressure  limit 
and  for  ES/RTS  >>1,  the  leading  terms  in  the  temperature  sensitivity  are 
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The  temperature  sensitivity  for  the  BDP  monopropellant  model  has  been  previously 
derived, 13,30  Although  it  has  not  been  expressed  in  the  same  form  as  Eqn.  (59) 
above,  the  two  expressions  can  be  shown  to  be  equivalent  except  that  the  third 
term  on  the  RHS  of  Eqn.  (61)  is  absent  in  the  BDP  formulation  because  of  a 
slight  difference  in  definition  of  x*  between  Ref.  11  and  Refs.  12,  13  and  30. 
Coates^l  has  obtained  an  expression  of  a  based  on  a  laminar-flame-theory 
model.  Its  dominant  term  is  equivalent  ?o  the  second  term  on  the  RHS  of  Eqn. 
(59).  In  laminar-flame-theory  models3>15, 26, 31  the  burning  rate  is  determined 
entirely  from  the  conditions  prevailing  in  the  thin  reactive-diffusive  zone 
at  x*.  There  being  no  explicit  dependence  of  M  on  T0  or  Ts,  only  the  implicit 
dependence  of  Tf  on  TQ  remains  to  govern  the  temperature  sensitivity.  This 
is  the  same  source  of  the  last  two  terms  in  Eqn.  (59).  As  previously  noted 
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the  dependence  of  these  terms  on  TQ  is  qualitatively  different  than  that  of 
the  empirical  Eqn.  (52). 

3.  QCHR  Approximation  to  <j>g 


Requiring  Eqns.  (18)  and  (47)  to  be  equal  and  utilizing  Eqn.  (19)  leads 
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in  the  high  pressure  limit  and  the  constant  value  a^  (Section  IVA1)  for 
sufficiently  low  pressures.  The  pressure  exponent  indicated  in  Eqn.  (60) 
is 
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(61) 


for  high  pressures.  It  is  of  considerable  interest  that  this  pressure  index 
should  be  explicitly  coupled  with  the  solid  and  gas  phase  activation  energies 
as  well  as  the  reaction  order.  An  expression  such  as  Eqn.  (61)  has  potential 
value  to  performance-tailoring  efforts.  It  also  exemplifies  the  value  of  an 
analytically  tractable  idealization  in  making  a  connection  between  practical 
burning  rate  behavior  and  more  fundamental  processes  which  might  be  characterized 
independently  of  burning  rate. 

The  variation  of  Ts  with  T0  can  be  shown  to  be 
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Since  (2Es~EG)  is  of  order  Es>  in  the  high  pressure  limit  and  for  »1, 
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C  is  close  to  unity.  Therefore  the  leading  term  in  a  is 
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Thus  the  temperature  sensitivity  is  also  seen  to  depend  explicity  on  the 
solid  and  gas  phase  activation  energies  and  in  the  same  mutual  relationship 
as  for  the  pressure  exponent.  Since  Ts  depends  strongly  on  Es  (which  can  be 
seen  by  equating  Eqns .  (19)  and  (60)  and  solving  for  Ts)  ,  the  QCHR  model 
suggests  that  one  might  be  able  to  lower  both  the  temperature  sensitivity  and 
the  pressure  exponent  of  a  propellant  by  increasing  its  Es  (providing  of  course 
that  the  gas  phase  reactions  retain  similar  activation  energies  and  the  surface 
heat  release  does  not  change  much) . 

B.  Evaporation  Decomposition  Mechanism 

In  adapting  the  evaporation  law  (Eqn.  20)  to  a  Level  II  treatment  (i.e., 
no  diffusion)  we  assume  that  mg-0  =  1  so  that  Pg  =  P.  Thus  Ts  is  a  function 
of  P  alone.  Also  Qs  must  be  reinterpreted  as  the  heat  of  transformation  (Qv) 
and  will  have  the  opposite  sign  as  T0,  as  for  example  in  the  redefinitions 
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Since  <f) s  now  depends  on  P  (through  Ts)  as  well  as  T0,  it  is  no  longer 
possible  to  effect  a  separation  of  the  influences  of  T0  and  P  as  was  done  in 
the  pyrolysis  case.  Furthermore,  making  the  reasonable  assumption  that  the 
surface  temperature  will  never  fall  below  T0,  it  is  evident  that  <f>g  does  not 
go  to  zero  at  low  pressures  (due  to  sign  change  on  Qv)  so  that  the  burning 
rate  will  not  approach  a  constant  at  low  pressures  as  was  found  for  the 
pyrolysis  mechanism.  The  burning  rates  and  temperature  sensitivities  for 
the  different  models  are  obtained  as  follows. 

1.  CTRR  Approximation  to  <$>q 

The  burning  rate  for  the  evaporation  mechanism  coupled  with  the  CTRR  <$>q 
is  obtained  by  equating  Eqns.  (18)  and  (53)  and  solving  for  M.  One  obtains 
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It  can  be  expected  that  Ev/RTs  will  normally  be  large  so  that  the  pressure 
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dependence  of  the  bracketed  quantity  will  be  dominated  by  pv^  over  the  entire 
pressure  range. 

Since  Ts  does  not  depend  on  T0,  Op  may  be  found  directly  by  applying 
Eqn.  (51)  to  Eqn.  (63).  The  result  is 
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2.  DFHR  Approximation  to  <f> 

For  an  evaporation  mechanism  the  expression  for  M  is  similar  to  Eqn.  (58) 
except  that  Ts  depends  directly  on  P  (through  Eqn.  20),  and  Qs  is  appropriately 
redefined.  The  temperature  sensitivity  is  obtained  directly  from  the  revised 
Eqn.  (58).  The  result  is  the  same  as  Eqn.  (59)  with  appropriate  redefinitions 
of  Ts,  Qs,  and  QT. 


3.  QCHR  Approximation  to 

VJ 

Again  equating  Eqns.  (18)  and  (47)  but  this  time  using  the  evaporation 
law  (Eqn.  20),  the  following  burning  rate  is  derived 


(65) 


appropriate  to  all  pressures.  The  temperature  sensitivity  for  this  case  is 
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C.  Discussion 


Several  comments  are  in  order  concerning  the  formulas  just  derived. 

First  of  all,  in  discussing  pressure  dependences  we  have  consistently  taken 
terms  of  the  order  of  Ts  as  slowly  varying  provided  ES/RTS  >>1.  This  was  done 
in  order  to  obtain  the  essence  of  the  functionalities  involved.  In  all  of 
the  high  pressure  expressions  for  M  and  ap,  Ts  appears  in  the  denominator. 
Since  by  either  the  pyrolytic  or  evaporative  mechanism  Ts  will  increase  with 
increasing  pressure,  M  will  tend  to  increase  with  P  somewhat  less  rapidly 
than  the  dominant  pressure  exponent  would  suggest.  Typically,  calculations 
in  which  Ts  is  obtained  by  solving  the  transcendental  equation  <|>g  =  for  a 
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pyrolysis  mechanism  resulted  in  a  pressure  exponent  about  10-20%  lower  than 
that  given  by  Eqn.  (65).  However,  the  Pn  dependence  of  M  is  preserved, 
i.e.,  n  is  constant  to  a  very  good  approximation  in  the  high  pressure  limit. 

Such  an  accuracy  is  well  within  the  range  of  usefulness  in  view  of  the  consider¬ 
able  uncertainties  involving  the  kinetic  activation  energies  and  possibly  the 
model  idealization  itself. 

We  should  also  remark  on  our  neglect  of  radiative  losses  from  the 
surface  of  the  solid.  This  radiative  loss  term  would  be  added  to  the  RHS 
of  Eqn.  (18)  and  would  alter  the  character  of  4>s  as  Ts  ->•  Ts '  at  low  pressures. 

At  values  of  P  sufficiently  low  that  the  radiative  loss  term  dominates  <J>  an 
abrupt  deflagration  limit  (below  which  no  steady-state  solution  exists)  may 
result  as  has  been  discussed  by  Johnson  and  Nachbar^  and  Buckmaster,  et  al . 
Crawford,  et  al.,21  notes  such  a  apparent  limit  at  about  1  atm  for  a  double¬ 
base  propellant  (offering  no  explanation  for  it).  Above  this  pressure  they 
report  that  the  burning  rate  is  well  described  by  Eqn.  (50)  with  a  non-zero 
value  of  a.  The  low  pressure  trend  toward  a  constant  value  of  M  is  therefore 
observed  but  the  limit  itself  is  probably  masked  by  the  eventual  dominance  of 
the  radiative  losses.  It  is  also  possible  that  the  apparent  decrease  of  Qs 
with  decreasing  pressure  observed  in  thermocouple  experiment sl7, 33  at  low 
pressure  may  be  due  to  radiative  losses  and  not  a  change  in  reaction  mechanism. 

Collecting  the  results  of  this  section  we  find  that  all  of  the  gas-phase 
models  considered  lead  to  a  burn  rate  of  the  form  M=a+bPn  for  a  pyrolysis 
mechanism  and  M=bPn  for  an  evaporation  mechanism.  It  is  at  least  possible, 
therefore,  to  account  for  the  experimentally  observed  functional  forms  of  M(P) 
using  only  a  single  gas-phase  reaction.  However,  derivation  of  these 
functional  forms  are  only  a  necessary  and  not  sufficient  condition  for 
verifying  the  idealization,  as  is  readily  apparent  from  the  diversity  of  gas 
phase  approximations  leading  to  the  same  functions  for  M(P) .  This  circumstance 
is  simply  a  formal  articulation  of  the  widely  recognized  fact  that  a  more 
detailed  experimental  constraint  than  the  burning  rate  pressure  dependence  is 
needed  to  provide  an  effective  test  of  a  proposed  combustion  model . 

The  temperature  sensitivity  functions  proved  more  responsive  to  the  gas 
phase  approximation  used.  Only  the  QCHR  approximation  led  to  a  form  of  <jp 
identical  with  the  empirical  form.  The  other  two  approximations  are  close, 
however,  differing  only  by  two  additional  terms  which  depend  inversely  on  Tf. 
These  extra  terms  are  in  general  probably  not  negligible  and  may  represent  a 
significant  variation  of  the  CTRR  and  DFHR  expressions  for  Op(T0)  from  the 
empirically  deduced  Eqn.  (52). 
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V.  SUMMARIZING  REMARKS 


The  steady  combustion  of  homogeneous  solid  propellant  is  likely  to  involve 
a  myriad  of  complexities  whose  details  depend  intimately  on  the  particular  pro¬ 
pellant  being  studied.  In  spite  of  this  inherent  complexity,  most  of  these 
propellants  burn  with  a  rate  whose  pressure  dependence  can  be  well  represented 
empirically  by  one  or  the  other  of  these  simple  forms 

r  =  a  +  bP 
r  =  bP™ 


Their  temperature  sensitivities  are  commonly  described  by  the  empirical  form 


a 

P 


c 

d-T 

o 


We  have  shown  that  these  same  formulas  are  explicitly  derivable  by  assuming 
that  the  burning  process  consists  of  a  single  condensed-phase  process  and  a 
single  gas  phase  reaction.  If  the  condensed-phase  process  is  assumed  to  be 
endothermic  (such  as  a  phase  change) ,  the  idealization  leads  to  the  form 
r=bPn.  For  an  exothermic  condensed -phase  process  (such  as  a  pyrolysis  re¬ 
action)  the  idealization  leads  to  the  more  general  form  r=a+bPn  (which  includes 
the  possibility  that  "a"  may  be  negligibly  small).  These  same  results  for  r(P) 
were  obtained  even  when  rather  different  approximations  to  the  heat-transport 
problem  in  the  gas-phase  were  used.  However,  only  for  the  new  gas-phase 
approximation  (QCHR)  developed  in  this  study  was  a  temperature  sensitivity 
formula  recovered  which  is  identical  to  the  most  common  empirical  form. 

It  was  shown  that  many  of  the  existing  gas-phase  approximations  could 
be  categorized  into  one  of  three  basic  assumptions  regarding  the  behavior  of 
the  heat  release  function.  This  enables  one  to  test  the  relative  suitability 
of  these  assumptions  for  a  given  propellant  by  comparison  with  empirical 
heat  release  profiles  calculated  from  temperature  profiles.  Such  an  approach 
has  the  advantage  that  no  a  pnLofvL  estimates  of  kinetics  constants  are  needed. 
The  comparison  was  made  using  published  embedded  thermocouple  records  for 
nitrate  ester  propellants.  Both  the  delta  function  heat  release  (DFHR) 
approximation  and  the  quasi-constant  heat  release  (QCHR)  approximation 
appeared  to  be  roughly  consistent  with  the  heat  feedback  one  would  compute 
from  the  transformed  temperature  profiles. 

With  regard  to  the  applicability  of  the  basic  idealization  considered, 
we  have  shown  only  that  it  is  capable  of  accounting  for  the  observed  shapes 
of  r(P)  and  0p(To).  The  significance  of  this  conclusion,  we  feel,  is  that 
it  was  reached  largely  without  recourse  to  specific  values  of  the  reaction 
rates,  such  rates  being  a  matter  of  considerable  controversy  at  the  present 
time.  On  the  other  hand,  the  derived  expressions  for  burning  rate  and 
temperature  sensitivity  differ  in  detail  among  the  various  models.  Without 
concrete  values  for  all  the  parameters  one  cannot  judge  the  extent  of  this 
disagreement.  For  this  reason  the  first  task  of  a  continued  effort  should  be 
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to  determine  the  relative  and  absolute  accuracy  of  these  models,  possibly  by 
comparison  with  numerical  calculations  for  a  few  plausible  data  sets.  In 
addition,  further  attention  might  profitably  be  given  to  assessing  and 
extending  the  reliability  of  the  embedded  thermocouple  technique  as  well  as 
exploring  techniques  for  extracting  further  information  regarding  the  kinetics 
of  the  fizz  zone  reactions  from  the  temperature  profiles. 
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B  -  gas-phase  reactant  label 

C  -  gas-phase  product  label 
B  C  C 

Cp>  C  ,  Cp  -  specific  heats  at  constant  pressure  in  gas-phase 
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u  -  local  mass  flow  velocity  in  gas  phase 
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g  -  term  in  -^=r-  expression  for  DFHR  model 

o  1 

0 
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Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports . 

1 .  BRL  Report  Number^ _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name: 


Telephone  Number: 
Organization  Address: 


